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The primordial spectra are calculated using dispersion relations which deviate from the relativistic 
one above a certain energy scale A. We determine the properties of the leading modifications with 
respect to the standard spectra when A ^ H, where H is the Hubble scale during inflation. To 
be generic, we parameterize the lowest order deviation from the relativistic law by a, the power 
of P/A where P is the proper momentum. When working in the asymptotic vacuum, the leading 
modification scales as {H/A)°' for all a, except for a discrete set where the power is higher. Moreover, 
this modification is robust against introducing higher order terms in the dispersion relation. We 
then algebraically deduce the modifications of scalar and tensor power spectra in slow roll inflation 
from modifications calculated in de Sitter space. The modifications do not exhibit oscillations unless 
the dispersion relation induces some non-adiabaticity near a given scale. Finally, we explore the 
much less studied regime where H and A are comparable. Our results indicate that the project of 
reconstructing the inflaton potential cannot be pursued without making some hypothesis about the 
dispersion relation of the fluctuation modes. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

In agreement with the observed temperature 
anisotropies in the Cosmic Microwave Background 
inflation predicts an almost scale invariant spectrum 
of primordial density perturbations 0, These adi- 
abatic perturbations arise from the amplification of 
vacuum fluctuations of linearized metric and inflaton 
perturbations. This mechanism relies on a semi-classical 
description, i.e., quantum fields propagating in a 
background spacetime. However in most models, the 
inflationary phase lasts so many e-foldings that the 
fluctuations we today observe stem from vacuum fluctu- 
ations characterized by wavelengths much shorter than 
the Planck length at the onset of inflation. In this regime 
the semiclassical description is no longer trustworthy. 
It is therefore of importance to find out to what extent 
the inflationary predictions depend on the physics that 
takes place above (or near) the Planck scale. 

In the absence of a theory of quantum gravity, there is 
no obvious way to address this question. In fact several 
approaches have been adopted. The flrst one is directly 
inspired by what was done in black hole physics where a 
similar problem occurs [1, d, @] : to test the robustness of 
the predictions, one introduces some dispersion above a 
certain energy scale A, and studies the sensitivity of the 
power spectrum to this modification 0, Q • It was shown 
that the properties of the power spectrum are robust, i.e., 
the deviations with respect to the standard spectrum are 
small whenever the adiabaticity of the evolution is pre- 
served and A taken well above the Hubble scale H during 
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infiation [9*]. However, the precise relationship between 
the modifications of the dispersion relation and the in- 
duced modifications of the spectrum was not obtained in 
the general case. To get a generic estimate of the modi- 
fications a simplified approach was proposed in Ref . [l3] ■ 
Instead of introducing dispersion around the scale A, the 
vacuum state was imposed when the momentum P = A 
rather than in the asymptotic regime P — > oo as done 
when working in the Bunch-Davies vacuum. A third ap- 
proach based on an effective action obtained by integrat- 
ing out heavy degrees of freedom of characteristic mass 
A, was also used in Ref. [11] and a different estimate of 
the modifications was obtained. 

If all approaches agree on the fact that the inffationary 
predictions are robust when A/H 3> 1, there is indeed no 
agreement about the general properties of the modifica- 
tions of the spectra. In particular, there is no agreement 
concerning the power of H/A which characterizes their 
amplitude. In Ref. [IQi] it was argued that the devia- 
tions should generically be first order in H/A, whereas 
in the effective lagrangian approach it was argued that 
the corrections should be (at least) second order and con- 
tain only even powers of H/A. Moreover in a reformula- 
tion of the model of [l^ where the state imposed when 
P = A is the (properly defined) instantaneous adiabatic 
vacuum T2], it was shown that the corrections are third 
order. As of the first method based on dispersion, as we 
said, no general result seems to have been obtained. 

In addition, besides the question of the amplitude of 
the deviations, there is no agreement either on the generic 
character of the rapid oscillations which have been found 
in m, [11, [31 (but not in [ll| ) , confronted with observa- 
tional data in [H, [l^ , and criticized in [13] ■ 

In the present work, we aim to settle the questions con- 
cerning the amplitude of the deviations and the presence 
of fast oscillations when using dispersion, and still work- 
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ing with the asymptotic (Bunch-Davies) vacuum. To get 
results which are not bound to a particular dispersion re- 
lation (or to a particular class thereof), we parameterize 
the first deviation with respect to the relativistic law by 
a scale A and a power a in the following way 

= F\P^) = P^(l± (P/A)" + O ((P/A)'")) . (1) 

The sign determines whether the propagation is super- 
(+) or subluminous (— ). As in former works, the pre- 
ferred frame which is used to define and P is taken 
to coincide with the cosmological frame: fl is thus the 
proper frequency and P the norm of the spatial momen- 
tum as measured by comoving observers. In these mod- 
els the isotropy and the homogeneity of FLRW are pre- 
served, but the local 0, [l^ Lorentz invariance has been 
broken. (This is not the case in [HI.) It should also 
be noticed that these models respect a modified (weak) 
Equivalence Principle fl^ in that, when considering high 
momenta (in the preferred frame) and neglecting the 
gradients of the metric, the physics is the same as in 
Minkowski space (in the preferred rest frame). 

To separate the modifications due to the dispersion 
above the scale A from those governed by slow-roll pa- 
rameters, we first compute the deviations of the spectrum 
of a scalar field propagating in de Sitter space and then 
show how these determine, by simple substitutions, those 
of scalar and tensor modes in slow-roll inflation. 

In Section |TT1 to start the analysis, we algebraically 
solve a particular case; a quartic superluminous disper- 
sion (a = 2 in eq. H])), in do Sitter space. We compute 
the resulting spectrum for all values of H/A. This ana- 
lytical treatment allows us to identify the nature of the 
signatures, and to prepare the numerical treatment. 

In Section [ml using numerical integration techniques, 
we consider the general case, with a ranging from 1 to 6. 
When H/A ^ 1, we first show that the dominant devi- 
ation of the power spectrum scales as (H/A)". In other 
words the leading modification is linear in the lowest or- 
der deviation of the dispersion relation. This is true for 
both super- and subluminous dispersion, and for all val- 
ues of a but for a discrete set of powers (a^ = 3-1- 2i) 
where an overall a-dependent factor vanishes and where 
the leading modification scales with a power higher than 
ai. Secondly, we show that higher order terms in the 
dispersion relation are irrelevant in that they induce sub- 
dominant deviations which vanish faster than (H/A)" 
in the limit H/A 0. The signatures are thus ro- 
bust against modifying the dispersion relation by adding 
higher order terms. Thirdly, we show that the signa- 
tures of super- and subluminous dispersion have equal 
magnitude and opposite sign for all a ^ at, and all val- 
ues of H/ A ^ 1 . Fourthly, when the dispersion relation 
is smooth and the asymptotic vacuum well-defined, the 
modifications of the spectrum do not display oscillations. 
On the contrary, (fast) oscillations appear when some 
non-adiabaticity is localized near a certain (UV) scale, 
in agreement with the conclusions of Ref . [17| . We also 



briefly study the other regime where H/A is comparable 
or greater than 1. (For subluminous dispersion we chose 
the asymptotic behavior of F so that the the adiabatic 
vacuum is well deflned.) As one might have expected, 
the deviations are no longer governed by the lowest or- 
der modification of the dispersion relation. This does 
not mean however that this regime should be discarded 
because it is not known at which scale the semi-classical 
description loses its validity. (In some higher dimensional 
models (see [l^l and references therein), this scale could 
be much smaller than the (4 dimensional) Planck mass, 
and therefore could be smaller than H.) 

We then establish in Section [TVl that the modifications 
of the spectra in slow-roll inflation can be obtained from 
the above results. The main result is that the modifica- 
tions depend on the wave vector k only through the ratio 
Hk/A, Hk being the value of the Hubble scale when the 
/c-mode exits the Hubble radius. 

In the Conclusions we briefiy discuss the possibility of 
distinguishing between modifications of the spectra com- 
ing from dispersion and other modifications, like those 
stemming from a change in the inflaton potential. In par- 
ticular, we point that our results predict a violation of 
the consistency relation, even in the case when the scalar 
and tensor modes obey the same dispersion relation. 

II. MODIFIED POWER SPECTRUM 

A. The model 

In this Section, we consider a minimally coupled mass- 
less scalar field $ propagating in a FLRW background 
spacetime in comoving coordinates: 

ds^ = a{T]f {-drf + d^) . (2) 

The spectrum of scalar perturbations and gravitational 
waves can be deduced from the power spectrum of this 
field, see d, [2l[ and Section [iVl We then introduce a 
non-linear dispersion relation that we parameterize by P, 
see eq. ([1]). We assume that one recovers the standard 
relation in the IR, i.e. p2(p2) ^ p2 -yy^g cafl A the UV 
scale which weighs the first non-linear term of P^. 

Decomposing the field in Fourier modes with fixed co- 
moving momentum k = aP, the equation of the rescaled 
mode = a^k reads 

Before considering slow-roll inflation, we flrst specialize 
to de Sitter space, where 




since a — —1/ {Hif) . All dependencies in k drop out when 
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using the variable x = —krj = P/H: 
dl 



(5) 



Thus, the power spectrum remains scale invariant in de 
Sitter space when the dispersion relation is expressed in 
terms of proper frequency and momentum, i.e., when it 
respects the Equivalence Principle, and when the pre- 
ferred frame coincides with the cosmological frame. [39[ 



B. Relativistic case 

Let us briefly recall the calculation of the power spec- 
trum for the standard relativistic case — . The 
mode equation ([5]) reduces to: 



C. Quartic dispersion relation 

In this subsection, we compute the power spectrum in 
the particular case 



Fi 



P' 1 



p2 

A2 



(10) 



The subscript 2 indicates that the first (and only) non- 
linearity in F is quadratic in P/A and the + sign indicates 
that the dispersion is super luminous. This dispersion 
relation 1241 has been studied in a cosmological context 
in Ref. [2^. However, to our knowledge, the following 
calculation has never been made. We believe this is the 
only dispersive case where an exact calculation of the 
power spectrum is possible in terms of hypergeometric 
functions. 



dl 



(6) 



To obtain the power spectrum one needs to identify the 
asymptotic positive frequency solution of the above equa- 
tion. Indeed, in any successful model of inflation, the 
relevant modes we today observe had proper momenta 
obeying P ^ at the onset of inflation, and hence were 
all in their ground state [l^]. The power spectrum P 
of the relevant fluctuations is thus given by the follow- 
ing VEV 

Piv,k) = ^ /'d3^e'*^«(0|$(r;,0d(r;,0)|0) 



27r2 



(7) 



where |0) is the Bunch-Davies vacuum [231 and where 
is the positive unit norm Fourier mode associated with 
this asymptotic state. 

The observationally relevant quantity is the value of 
P(fc) at late times, long after horizon exit, k/aH = x ^ 
0. When ^^^{rj) is written in terms of 



u'"(a;) 



1 
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the unit wronskian solution of eq. 
itive frequency for krj = —x ~c 



(8) 



3) that is purely pos- 
one obtains 



27r2Po = H- 



IP 

2 



(9) 



The index stands for the unperturbed relativistic case. 

Had we worked in slow-roll inflation rather than de 
Sitter space, Po(fc) would have been given by the above 
with H replaced by Hk , the value of H when the fc-mode 
exited the Hubble radius, i.e., k/{a}~H}~) — 1, see Sec- 
tion [IV] for details. 



Analytical expression of the power spectrum 



When F is given by eq. H 
comes: 



X 

4A2 



the wave equation be- 



(11) 



From the comparison of this equation with eq. ([6]), one 
can deduce that the modifications of observables with 
respect to the relativistic case will all be governed by the 
dimensionless parameter 



A 



A 
2H 



(12) 



(The factor 1/2 has been introduced for convenience and 
will be retained throughout the paper.) 

As in the relativistic case, the state of the field is cho- 
sen to be the Bunch-Davies vacuum, that is, the ground 
state associated with the asymptotic positive frequency 
solution of eq. pT|) . This solution, normalized to unit 
wronskian, is given by 



^2+ 



(x) 



(13) 



The definition of the Whittaker function W^,^ can be 
found in [2^, p. 505. That the function u^2+ u.nit 
wronskian and is positive frequency at early times can be 
seen from its asymptotic form for large x (see equation 
13.5.2 in [1^): 



^2+ 



(x) 



X 



+41 



n fx) 



(14) 



This asymptotic form is precisely the WKB solution of 
eq. (jlip with positive frequency. One verifies that the 
corrections vanish as l/x^ when a; — s- oo. Therefore the 
Bunch-Davies vacuum is well defined. 



The power spectrum is then straightforwardly obtained 
from the behavior of W^x 3 for a; ^ (equation 13.5.6 
in dl). This gives 



'5 _ ,-A^ 



X (l + 0(a;2)) . 
Using this expression we obtain 

27r2p2+(A) = H^(x'' 1 4"+ (a;) 



(15) 



7r(2A) 
X 



3/2, 
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(16) 



This expression is exact and vaUd for aU values of H/ K. 
It gives the power spectrum when using the dispersion 
relation of eq. (|10p and when working in the Bunch- 
Davies vacuum. 

Using equation 6.1.45 of [26], one verifies that the fac- 
tor of the relativistic spectrum Pq = i?^/2 tends to 1 
when A = K/2H -^00. In this we recover that the 
spectrum is robust, that is, the standard value obtains 
when adiabaticity is preserved, as it is here the case when 
i//A< 1. 




Figure 1: Plot of the power spectrum, eq. (|16|l . for the quartic 
superluminous dispersion relation eq. (|10p divided by Po, as 
a function of = 2H/A. On the left side, for A — + 00 
(A^^ 0), the relativistic spectrum is recovered, while for 
A ^ 0, the power is suppressed as A^^^. 



large A, in the adiabatic regime, the modified spectrum 
asymptotes to the standard result. Instead, for small A, 
the modified power spectrum tends to zero as 



P2+(A)~PoX 



47r 



(2A) 



3/2 



(18) 



2. Signatures for small H/ A 

Since we have the exact expression of the power spec- 
trum, we can properly extract the first corrections in 
H/A = 1/2A. When A > 1, using equation 6.1.47 of [H], 
the r.h.s. of (|16p can be expanded into: 



P2+(A) =Po 1 



+ 0(A-4) (1 + e--^) . (17) 



This expression contains three factors: the relativis- 
tic spectrum and two factors which tend to one when 
A ^ 00. The first one contains a polynomial starting 
with a quadradic correction in 1/A, whereas the correc- 
tion term of the second is exponentially suppressed. In 
Appendix |^ we show that an approximate calculation 
based on WKB waves correctly reproduces these features, 
namely that there exist two sources of corrections: poly- 
nomial corrections related to local modifications of the 
mode with respect to the relativistic case, and exponen- 
tially small corrections related to non-adiabatic transi- 
tions. 



In Figure [21 to display the signatures of the disper- 
sion, we analyze the relative deviation with respect to 
the relativistic spectrum, AP2+/P0 = (P2+ — Po)/Po. 

In the upper plot we have represented log|AP2+/Po| 
as a function of logA~^. One clearly sees that when 
1/A = 2H/A < 0.1, the polynomial correction term 
in pT]) largely dominates all other contributions. In- 
deed the curve asymptotes to a straight line, with the 
expected slope, equal to 2, and the expected intercept 
= log 5/16 ~ —0.505. In the lower plot, in anticipation 
to the numerical analysis, we have represented the quan- 
tity 



Af(A) = A' 



Pf-Pq 
Po 



(19) 



where a is the power of P/A in the first non-linearity in 
F, a = 2 here. The interest of this representation is that 
as soon as the deviation in A~" becomes the dominant 
one, the curve becomes horizontal. 



III. NUMERICAL ANALYSIS 



3. Figures 



A. Parameterization of the dispersion relation 



To visualize the signatures and to prepare the numeri- 
cal analysis of the next Section, we have represented sev- 
eral figures. In Figure [1] we show P2_|_ of eq. ^TE\i divided 
by Po as a function of A"^ for 10^^ < A^^ < 10^. For 



In the previous section, when studying the dispersion 
relation with a quadratic correction term in P/A, we 
found that the first deviation of the power spectrum 
was quadratic in H/A. Our first aim is thus to verify 
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whether the power of H/ K characterizing the first devi- 
ation is always the same as that of P/K in the leading 
non-linear term of the dispersion relation. Secondly, we 
want to show that when adding subleading correction 
terms to the dispersion relation, i.e., terms character- 
ized by a power of P/A higher than a, the signatures 
are robust, i.e., they are insensitive to these additional 
terms. Third, we want to analyse in a symmetrical man- 
ner super- and subluminous dispersion. Therefore, we 
select subluminous dispersion relations which possess a 
well defined asymptotic regime P — )■ oo so that the adi- 
abatic vacuum condition can be implemented for early 
times X = —ki] oo. (Remember that subluminous 
dispersion relations containing only one non-linear term 
do not possess such a well defined regime, because ft^ 
becomes negative for P/A > 1.) 

To these ends, we consider the following parameteriza- 
tion of the dispersion relation: 



^2 ^ F'^{P,A;a,P,N) 
= < /3P + 7A tanh ° 



(20) 



CP 
A 



p n2q + 2 



where the extra parameters C, 7 verify the following equa- 
tions: 



4a|l-/3| 



l/a 



1-/3 



(21) 




log X 



Figure 2: Plot of the relative corrections AP2+/P0, obtained 
from the exact expression (|16|l . In the log- log representa- 
tion (upper plot), the dominance of the first signature is seen 
through the linearity of the curve: when < 0.1 the curve 
is indistinguishable from the straight line with slope 2 and 
intercept log(5/16) ~ —0.505 (dashed line) associated with 
the quadratic correction in eq. (|17p . In the lower plot, it 
is seen through the constantness of A2+ — A^AP2+/Po for 
< 0.1, with a value equal to the coefficient of the first 
signature term in (fT7)l. —5/16 = —0.3125. 




When these are verified, the Taylor expansion of P^ for 
P/A < 1 is 



(22) 



We obtain a superluminous relation when (3 > 1 and a 
subluminous one when /? < 1. In addition, in each sector, 
at fixed a, varying /? modifies only the coefficient of the 
subdominant deviations (with powers of P/A equal to or 
greater than 2a). 

Returning to the exact expression (PO)) , we notice that 
dplnF ^ in the limit P/A 00. Therefore, in all 
cases we reach an adiabatic regime which allows to prop- 
erly define the asymptotic vacuum, see below for more 
details. The exponential factor in eq. (|20p might seem a 
priori useless since it does not affect any of the first three 
terms in eq. (|22|) . It has been introduced to facilitate the 
numerical integration because it improves the validity of 
the WKB approximation for large x. The initial vacuum 
conditions can then be specified for smaller values of x 
(typically .t;,! ~ 500A/C), thus avoiding a too large accu- 
mulation of numerical error coming from the (physically 
irrelevant) early evolution where P/A ^1. In the nu- 
merical integration we shall put N — 100/C — Xin/5A, 
thereby switching on the tanh about two e-foldings af- 
ter having started the integration. (We have checked the 
stability of the results when using higher values of N and 
Xi„/A.) 

Figure [3] represents the dispersion relation ((20)) for 
a — 2 and several choices of /3. In the left plot the 
cases P = 0.2 and f3 — 1.8 illustrate the fact that the 
parameterization of eq. (j20p encompasses both super- 
and subluminous dispersion relations. In the right plot, 
the subluminous dispersion relations with f3 — 0.1 and 
f3 — 0.5 are represented along with the quartic disper- 
sion relation = P^(l — P^/A^). Since the leading 
(quadratic) deviation is the same for all three cases, the 
curves coincide even after having left the linear regime, 
until P ~ A/2. 

To study the adiabaticity, we write the corresponding 
wave equation in the variable x = —kr] as 



(23) 



with the effective square frequency 
1 



F^iH'x',A;a,(3,N) - — . (24) 



The evolution is adiabatic whenever ap given by 

. (9a; In 



apix) 



UJF 



(25) 



(a) 



(b) 





(P/A) 



Figure 3: Plots of the dispersion relation eq. (I20|l for q = 2 and various values of f3. The linear relativistic dispersion relation 
is also represented in both plots, as well as the quartic subluminous dispersion relation in the right plot. In plot (b), Q,^ is 
shown as a function of so that the quartic dispersion relation be simply a parabola. 



remains much smaller than 1 



B. Numerical resolution 



In the absence of dispersion, the non-adiabaticity only 
arises from the — term and ctq decreases for x ^ 1 as 
When adding dispersion, the non-hnear behavior 
of F brings another source of non-adiabaticity. However, 
when x/A 3> 1 is also satisfied, the non-linearities intro- 
duced by the hyperbolic tangent in (|20[) are suppressed by 
the exponential factor. The remaining terms in ujf are 
a constant term (= /3^) and the (small) negative term 
—Ijx^. Thus in this asymptotic regime, op decreases 
as [ji'^x^)~^ thereby guarantying that the adiabatic vac- 
uum stays as well-defined as in the absence of dispersion. 
(Had we not introduced the exponential factor, the hy- 
perbolic tangent would asymptote to a constant which 
would bring a contribution \o ap decreasing only as l/x^, 
thereby imposing to use a larger value of to describe 
with a high precision the asymptotic vacuum.) 



The behavior of ctf is represented in Figure H] for sev- 
eral values of the parameter A and for subluminous (left 
plot) and superluminous (right plot) dispersion. The evo- 
lution is globally more adiabatic in the superluminous 
case, as expected since in this case is greater than 
in the subluminous case. For each value of A, there is a 
local maximum around x ~ A, whose height diminishes 
when A increases. When scale separation is not realized 
(A — 1), in the superluminous case a — 2 and (i — 1.8, 
this local maximum merges with the rapid increase oi ap 
when X approaches 1. In the subluminous case a = 2 
and /? = 0.2, dp takes non-negligible values, of order 
10""'^ around a; ~ 6, significantly before horizon exit. 



Following conventional techniques, the second order 
differential equation (|23p is separated into a system of 
4 first order equations, 2 for the real part of and its 
derivative, and 2 for the imaginary part of and its 
derivative. After setting the initial conditions, discussed 
below before eq. , this system is integrated using the 
embedded 8th order Runge-Kutta-Prince-Dormand algo- 
rithm provided in the GNU Scientific Library, from Xi to 
some Xf long after the Hubble scale exit (x = 1), typi- 
cally x/ — 10"'^. In this algorithm, the error is estimated 
at each step and the stepsize adapted to keep the errors 
within fixed bounds. The use of this algorithm was nec- 
essary, because the relative deviation APp/Pg reaches 
very small values when A is large. For instance it is of 
order 10~^ for a = 4 and A — 100. Since the high values 
of A also require a larger number n of integration points 
(because Xi ~ 500A/C, see below, is further away), the 
absolute error for each step must be kept under \^~^^ jn 
if we want to be able to distinguish values of APp / Pg as 
low as lO^i*'. 

The initial conditions are fixed using the fact that the 
WKB solution of eq. (|23p. see eq. (jA2p . becomes ex- 
act when x — > cxD. For some finite Xi ^ A, the differ- 
ence between the WKB solution and the exact positive 
frequency solution is of the order of ctj;- [I4I. Thus the 
WKB approximation becomes excellent whenever both 
conditions x/A ^ 1 and a; ^ 1 are satisfied, since in this 
case ap ~ {0^x\^~'^. In practice, the initial conditions 
are fixed at Xi = 500 A for A > 10 and at Xi — 5000 for 
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Figure 4: Influence of A on the adiabaticity of the evolution. 

ap is given in eq. (|25|l . In all plots a = 2. To get comparable values we have represented y/Xap- ctq for the relativistic 
dispersion relation is also shown in both plots for a comparison with the case A = 1. (a) Subluminous dispersion (/3 = 0.2). (b) 
Superluminous dispersion (/3 = 1.8). Subluminous dispersion is less adiabatic than superluminous, as might be expected since 
uiF is reduced in the first case. The local maximum at a; ~ A is higher for smaller values of A. 



A < 10, so that (Tp < 10 ^"7/3^. We thus safely impose 

1 



i UJp 



l + i 



1 dj:UJF 



(26) 



(The arbitrary phase of the mode has been put to at 
Xi without affecting the results.) 

At the end of the integration of the wave equation, the 
relative deviation of the power spectrum wrt the rela- 
tivistic spectrum is evaluated through: 



APf(x/) _ PF{xf)-Po 



= 2x}\4^^ixf)\' 



1 . 



(27) 



which follows from eq. ([9|). Since 2^/ is not exactly zero, 
the value obtained contains a finite contribution from the 
decaying mode. For A < 1, this contribution is negligible 
wrt the modifications due to dispersion and we simply 
ignore it. However, since for A 1 it decreases as x^ 
and since Xf ^ 10~^, the decaying mode contribution 
in (I27p is of order 10~^. It largely dominates the small 
corrections we are after. One could consider smaller val- 
ues of x/ but this leads to an increase of the number of 
integration points and of the numerical error. We thus 
use the fact that for large A and a; ^ 1, APf{x) is of the 
form: 



APFix) = Po{Af + Bfx''), 



(28) 



where Af and Bf depend on the values of the various 
parameters. Af is the contribution of the growing mode 



that we seek to extract. To do so, we compute the relative 
deviation ([27|) at 10a;/ and Xf. Using twice (|28|) . we 
properly extract Af- 



C. Results 

1. Properties of the signatures when H/A <C 1 

The expansion at large A of the analytical result of 
Section FlI C 11 and the qualitative arguments in Appendix 
1X1 led us to conjecture that APf/Pq is linearly related 
to the first deviation in the dispersion relation and 
thus possesses the following asymptotic form: 



APf 



6^ (a) 
A" 



for A 



(29) 



Since only the sign of /3 — 1 appears in the first devia- 
tion, is expected to depend on the superluminous (-1- 
exponent) or subluminous (— exponent) character of the 
dispersion, but not on the actual value of (3. 

In Figure [HI log APf / Pq is represented as a function 
of log A^^ for a series of subluminous dispersion relations 
with P — 0.2 and values of a from 2 to 4. For high values 
of A up to 10'^, we obtain a set of straight lines, with a 
slope precisely equal to a, thereby confirming the valid- 
ity of the asymptotic behavior given in eq. l(29)) . Had 
we continued the curves toward A = 1, we would have 
seen deviations from this linear behavior. This (non- 
adiabatic) regime around and beyond scale crossing will 
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Figure 5: Logarithm of the signature of dispersion for various 
powers a of P/A. (3 is fixed to 0.2. In each case, the slope is 
equal to the power a, see eq. (I22p . 



be studied in Section IIII Dl For superluminous disper- 
sion we obtain similar results which confirm the validity 
ofeq. ((29|) . 



0.5 



0.4 



0.3 



-0.3 



T ' 1 ' 1 ' r 
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|3=0.99 
|3=0.999 



P=1.001 
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log X 



Figure 6: Robustness of the signatures to higher order terms. 
a is fixed to 2. One verifies that, both for sub- (upper graph) 
and superluminous dispersion (lower graph), the deviations of 
the spectrum induced by the second non-linear term (which 
is weighted by (/3 — 1)~^) disappear for sufficiently large A to 
give rise to the leading deviation given by eq. 



2. Robustness of the signatures for H/A <^ 1 

To establish that higher order terms in P/A in the 
dispersion relation are irrelevant in the limit H/A ^ 1, 
in Figure [H] we have plotted Ai;- of eq. for a = 2 

and various values of /3 which weigh the term scaling 
as (P/A)^" in eq. (|22p . The value of /3 ranges from 
0.1, in which case the coefficient of this term is of or- 
der 1, to /3= 1 — 10^^ in which case the coefficient is 
~ 10^. Similarly for superluminous dispersion relations, 
f3 ranges from 1 -I- 10"'^ to 1.9. As expected, when A ^ 1, 
all curves agree irrespectively of the value of /3. We have 
numerically verified that the next order deviations {i.e. 
the departure from the asymptotic behavior of eq. ((29|) 
which are clearly seen in Figure scale as the square 
of {H/A)". We have also verified that their normalisa- 
tion coefficient contains a contribution which is linear in 
(1 — and which comes from the second correction 

term in eq. (|22p . as well as a contribution which must 
arise from the square of the first correction term in that 
equation. 

Since the function A^^^ is asymptotically constant when 
A/if ^ 00, the asymptotic properties of the signatures of 
UV-dispersion are governed by the two functions i5^(a) 
of eq. ([29]) . We have verified that similar results hold for 
all a ranging from 1 to 6, except for a = 3 and a = 5, 
which are particular cases as explained below. 

3. Properties of 5{a) 

To further investigate the asymptotic behavior of the 
signatures, we have represented in Figure [7] the function 



Ap{X) of eq. (fTO|) for various values of a and for both 
sub- (/3 — 0.2) and superluminous (/3 — 1.8) dispersion 
relations. 

Besides confirming the fact that A^? is indeed constant 
when A — > GO, we learn from Figure [7] two rather unex- 
pected results. First, when comparing super- and sublu- 
minous cases for each value of a, we see that S^{a) — 
—S~{a). More precisely, we have numerically checked 
that the sum Aa.+ {X) + Aa.-{X) scales as A~" at large 
A, and thus vanishes when A — s- cx). This result allows us 
to consider the unique function 

(5(a) =sign(l-/?)(5=^(a) (30) 

which does not depend on (3. We also point out that we 
have not been able to find any analytical explanation for 
this simple fact. See Appendix \X\ for an attempt in this 
sense. 

Second, we see that the sign of S changes between 2.5 
and 3.5. In fact it flips sign exactly for a — 3 when the 
dimensionality of spatial sections is 3. To establish this 
second fact, we consider the generalization of eq. ((24)) for 
a de Sitter spacetime of spatial dimension d: 

(31) 

To obtain S{a,d), we constructed a large number of 
curves Aa.d{X) and for each curve found the asymptotic 
constant value. We repeated this for different values of 
d (including unphysical non- integer values). The result 
is shown in Figure [S] (We verified that, for every d, 
5~^{a,d) — — (5~(a, d), thereby extending the validity of 
this peculiar correspondance.) 

We find that 6{a, d) vanishes precisely when a is equal 
to d. We also find that after the first zero, it flips sign 
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Figure 7: Comparison of super- and subluminous signatures. 
The legend applies to both plots. In (a) the dispersion is su- 
perluminous with P — 1.8. In (b) the dispersion is sublumi- 
nous with /3 — 0.2. In (a) we represented — Af as a function 
of log A so that the asymptotic values reached in each plot can 
be compared at the center of the figure. We clearly see that 
S+{a) = -5-{a). 



again precisely dA, a — d + 2 (for d = 2 we even veri- 
fied that it also vanishes at a = 6 = o? + 4) . We there- 
fore conjecture that S{a,d) changes sign at the points 
ai = d + 2i, I e N. In the particular case of the disper- 
sion relation PH)) . where a — 2, the exact calculation of 
Section III CI can be generalized to d spatial dimensions 
and one finds for large A: 



2+,d 



Po,d 



dd^ 

3 



4 1 
1^ A2 



(32) 



Thus in this case one verifies analytically that the coeffi- 
cient d^{2,d) vanishes for d — 2. 
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Figure 8: Normalisation (5(a,d) of the first signature as a 
function of a, for several values of the spatial dimension d. 
For each d. ?){n.f]\ vanishes a.nd chariEres sipn at the nnints 
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Figure 9: Logarithm of the relative corrections to the power 
spectrum for a = 3 and d = ?>. At large A./H, all curves are 
linear with the same slope, equal to 6. The dotted line shows 
the straight line of slope 6 coinciding with the curve j3 = Q.l 
at large A/H. 



4- Particular case a = 3 when d = 3 

In this subsection we return to 3 dimensions and con- 
sider the particular case a — 3 for which S vanishes. In 
Figure [9] we have represented log APp /Pq as a function 
of log A^^ for a = 3 and several values of /?. 

From the linear behavior for small A~^, we see that the 
corrections are still given by a power law. However, the 
slope is now 6, thereby establishing that the signature is 
given by the square of (H/A)". This dominant power of 
A~-^, i.e., the slope, is stable against changes of f3, but its 
coefficient, related to the intercept, is not. We verified 
that this coefficient depends linearly on (1 — /3)~^, like 
the coefficient N of the second correction term in the 
dispersion relation (see eq. ([^ ). but that it is not simply 
proportional to N. Thus when d vanishes, the dominant 
deviation of the spectrum comes both from the square 



of the first non-linear term in the dispersion relation and 
from a linear contribution of the second non-linear term. 



5. Fast oscillations 

In various works [l^, [H, it was found that the 
deviations from the standard spectrum display oscilla- 
tions with a high frequency proportional to A/H. In 
these works the UV scale A was introduced through the 
choice of the initial vacuum which was imposed when 
^in = fc/ain = A (i.e. Xin — 2X in our notations). De- 
pending on the adiabaticity of the choice of the vacuum, 
the power in H/A of the norm of the deviations runs from 
1 [Tq to 3 [H, but in all cases, the deviations oscillate 



with a high frequency. In Ref. the origin of these 
oscillations was attributed to the instantaneous charac- 
terization of the initial state. It was also shown that the 
oscillations are suppressed when smearing the UV scale 
(completely or partially depending on the width of A be- 
ing larger or smaller than H). 

Given this, it is interesting to look for oscillations when 
the UV scale A is introduced through a dispersion rela- 
tion. In agreement with the conclusions of Ref. , we 
find that the deviations display no oscillations when the 
dispersion relation is smooth, and when the initial state 
is sufficiently close to the asymptotic adiabatic vacuum. 
When instead the dispersion relation possesses a kink or 
a bump near a given UV scale, fast oscillations appear 
in the deviations, [i^l We have verified this in several 
examples. 

From this we get two conclusions. On the one hand, 
fast oscillations of the type found in Refs. fldl. fl^. [isL KE^ 
can also be obtained with dispersion even when work- 
ing in the asymptotic adiabatic vacuum. On the other 
hand however, to obtain them requires some odd fea- 
ture in the dispersion relation well localized near a UV 
scale which will cause a non-adiabatic evolution near that 
scale. Thus fast oscillations (with significant amplitude 
when compared to the deviations) cannot be considered 
as a generic property of the deviations of the spectra en- 
gendered by modifying the physics in the UV sector. 



D. Power spectrum for H/A > 1 

In our parameterization of the dispersion relation the 
asymptotic vacuum is well defined for all values of H/A. 
Hence we can investigate the properties of the power 
spectrum in the much less often studied regime where H 
is comparable and greater than A. (To our knowledge, 
this was never done before with dispersion, see (20l . [27| 
for an analysis in the presence of dissipative effects.) It 
is also of value to note that the simplest approach to the 
trans-Planckian question where there is no dispersion but 
where the quantum state of the field is fixed at some scale 
A (see [Mill P), is not defined in the regime H/A > 1. 

Figure [TOl represents APp/Po as a function of log 
for one super- and three subluminous dispersion relations 
with a = 2. The main plot shows that the different spec- 
tra, that merge into each other for A > 10, depart from 
each other and depend strongly on A when A approaches 
1. However, for A ^ 1, all curves asymptote to a constant 
value that strongly depends on (3. As already pointed out 
in Section FlH Al the asymptotic constant behavior can be 
understood from eq. (|20|) : when A ^ 1 there is actually 
no dispersion before and at horizon exit, but the mod- 
ified speed of the modes is /? ^ 1. Thus, when deahng 
with eq. (PU)) . the power spectrum does not depend on 
A and is equal to Pq / Z?^ (for all f3, lower or greater than 
1), in agreement with what we see in Figure [TUl 

Thus, as expected, when the condition A ^ 1 is not 
satisfied, the knowledge of the full dispersion relation is 
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Figure 10: Corrections to the power spectrum when scale 
separation is not realized, a is fixed to 2. The insert shows 
that the deviations for /3 = 0.1 and /3 = 0.2 are very different 
for A"^ > 1 even though they agree for A"^ < 1. 



needed to make predictions. However, the power spec- 
trum is still well defined as a VEV when the asymptotic 
vacuum makes sense. Therefore there is a priori no rea- 
son to exclude the possibility that the regime H > Ahe 
relevant for observational cosmology. And were this the 
case, this would complicate the identification of the slow 
roll parameters, not to mention the reconstruction of the 
inflaton potential. 



IV. 



SLOW-ROLL INFLATION 



To show that the above results apply mutatis mutandis 
to slow-roll inflation, we need to say a few words about 
the kinematics of slow-roll and power law inflation. 

A large class of realistic quasi-de Sitter backgrounds 
can be described in the framework of the slow-roll ap- 
proximation where the evolution of the Hubble scale is 
characterized by the smallness of the parameters 



ei 
£2 



1 

H~ 



dAnH < 1 



2H 



9t In ei < 1 , 



(33) 
(34) 



The linear order slow-roll approximation consists in keep- 
ing only the terms that are linear in ei and £2, as well 
as treating these parameters as constants. This requires 
that the condition dtt2 = 0(e^, £2, £162) be fulfilled. In 
what follows, expressions in slow-roll inflation should be 
understood as given to linear order in ei and €2- For 
the particular case of power-law inflation, ei is constant, 
£2 — 0, and the scale factor is exactly given by: 



1 



1/(1-^1) 



(35) 



When £2 ^ 0, to linear order in the slow-roll parameters. 
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asi only depends on ei: 



asi 



-Hov 



l+ei 



(36) 



which agrees to (1351) to first order in ei. This equation 
follows from 



da 



1 + e 



aH 



'- + 0{eleleie2). (37) 



In slow-roll inflation, the time-dependent term in the 
square frequency of the tensor modes is —d^a/a, as in 
eq. ([3]), whereas that of scalar modes is —d^f/f with 
/ = o-V^i 0]- In power-law and in first order slow-roll 
inflation, these terms are always of the form — /i/?7^ where 
/i is a constant. In power-law inflation, the value of the 
constant is the same for scalar (S) and tensor (T) modes 
and is given by 



ei 



(l-^i)^ 



(38) 



For slow-roll inflation, the values differ and are respec- 
tively 



^J■T 



3ei , 
3(ei- 



(39) 
(40) 



Given the above discussion, we can deduce the prop- 
erties of dispersion-induced modifications of the spectra 
for scalar or tensor modes in slow-roll inflation, by simply 
comparing the mode equation to that of the test field in 
a de Sitter background. Moreover, when considering the 
leading modification of the spectra, since it is governed 
by the lowest order non-linear term of i^, we can restrict 
our comparison of the mode equations to the following. 
In de Sitter space, in d spatial dimensions, the relevant 
terms for a test field are 



dt 



i± 



(1;)° 



1 



0. 



(41) 



In slow- roll inflation, in 3 dimensions, using again the 
variable x = —krj, they are given by 



0fe = 0, 



(42) 



where the constant ^ is ^.pi, fir or fis- In power-law 
inflation, the new parameters are 



1-ei 



where Hh is 



H, 



api iVk 



Hp 
1 - e 



1/(1- 



(43) 
(44) 

(45) 



In first order slow-roll inflation, Asi(ei, k) is given by the 
linearization in ei of eq. 

Comparing equation (|42| with eq. (I4ip . one sees that 
the results of the previous sections apply with A replaced 
by A(ei, k), a hy ai, and d given by 



(46) 



The condensed notation d^ indicates a dependence on 
both ei and £2, and ^ can be /Zp/, ijt or us according 
to the case considered. Thus, if we denote by Po,e and 
Pa^e the relativistic and modified spectrum respectively, 
we have, in place of eq. ((29|) . 



S^{ai,de) 



when A(ei, k) 



Po,. A(ei,A:)"i ' 

oo. Using eqs. and (HJ), we get 



(47) 



-Pg.e ~ Pq.i 



S^iai,d,) 
(2-2ei)-"i 



Hk 

X 



(48) 



where the power of Hk/ K is a, independently of ei, £2. 

Eq. (|48p is the main result of the paper. For power-law 
inflation, it gives the leading modification to the power 
spectra to all orders in ei. For slow-roll inflation, it gives 
the leading modification at linear order in the slow-roll 
parameters. It applies both to scalar and tensor pertur- 
bations, with different expressions for de, see above. It 
establishes that the fc-dependence of the modifications to 
the spectra only arises through {Hk/A)". [4l| 

Given the above substitutions, the overall coefficient 
of (i/fc/A)" vanishes for the discrete set of values 



Using 



ai^ei = de + 2i . 



and (|40p. one finds 



ai^T.si = 3 - £1 + 2i , 
ai,s,si = 3 - £1 - e2 + 



2z, 



2i. 



(49) 



(50) 

(51) 
(52) 



For these values, the leading modification scales with a 
power higher than a^. 

A few remarks are in order here. Firstly, as long as 
Hk/A <C 1, one can extend the correspondance between 
de Sitter space and slow-roll inflation to an arbitrary 
number of terms of a general dispersion relation. In- 
deed, writing the n*'' term of the dispersion relation as 
c„(P/A)"", this yields c„(x/2A)"" in the wave equation 
in de Sitter space. Thus, as above, the corresponding 
term in slow-roll inflation will be given by A — *■ A(£i, fc), 
a„ — > an/(l — £1) (or q;„ ««(! + ei) in slow roll infla- 
tion). This correspondance term by term can be used to 
predict the power of A(£i, fc) for an arbitrary number of 
subleading modifications to the spectra. 

Secondly, in the particular cases where the dispersion 
relation is of the type ([^0)1 when N 00, i.e. where 



12 



all powers of P/A are multiples of a, the complete mode 
equation for slow-roll inflation can be obtained by making 
the above replacements directly in the effective frequency 
of eq. Thus in this case, the complete modification 

to the spectra in slow-roll inflation can be obtained from 
the de Sitter result, for arbitrary values of Hk/A. 



V. CONCLUSIONS 

In this paper, we have determined the properties of 
the modifications of the inflationary power spectra en- 
gendered by dispersion when the quantum state is the 
asymptotic vacuum. When the leading non-linear term 
in the dispersion relation is (P/A)", see eq. ([2^ . we 
have established that the following holds in the regime 
H/A < 1, 

• The leading modification of the spectrum behaves 
as S{a){H/A)", for all values of a, except for a 
discrete set of values where d vanishes. 

• This leading modification is robust in the sense that 
when adding to the dispersion relation terms con- 
taining higher powers of P/A than a, these give rise 
to subdominant modifications of the power spec- 
trum, as clearly shown in Figure [S) 

• For all values of the power a, except those where 
S vanishes, the leading modification of the spec- 
trum has equal magnitude and opposite sign when 
comparing sub- and superluminous dispersion, see 
Figures [7] and m 

• When the state is the asymptotic vacuum, the mod- 
ifications of the spectra display no oscillations, un- 
less the dispersion relation has some localized bump 
that would engender some non-adiabaticity at that 
scale. 

• These results, which have been derived in de Sitter 
space, apply to slow-roll and power-law inflation, 
both for scalar and tensor modes, by replacing H 
by Hk — k/ak and adjusting some constants which 
appear in the mode equation, see eq. (I48|) . 

In brief, when Hk/A <C 1, when the evolution stays adi- 
abatic, and when the state is the asymptotic vacuum, we 
see that it is not necessary to exactly know the disper- 
sion relation, since only the first non-linear term matters. 
In particular, in slow-roll inflation, this implies that the 
dependence on k of the modifications only arises through 
{Hk/A)". We also learn that, when a is unknown, no 
general prediction concerning the properties of the sig- 
natures of high energy dispersion can be drawn even in 
the regime H/A <C 1. 

We have also analyzed the 'non-adiabatic' regime 
where H > A that might be relevant in some theories 
with large extra-dimensions. The modifications of the 



spectrum are still well defined and governed by the dis- 
persion (when working in the asymptotic vacuum) even 
though they are no longer governed by the lowest order 
non-linear term of the dispersion relation, see Figure 1101 

To conclude we discuss to what extent the modifica- 
tions to the power spectrum wc obtain could be dis- 
tinguished, given some observations, from some other 
change in the scenario, such as, for instance, a change 
in the inflaton potential. 

The main point is that the absence of (rapid) oscil- 
lations in eq. (|48p greatly complicates this identifica- 
tion, since smooth deformations of the spectra can be 
accounted for by various means. As a consequence, the 
project of reconstructing 'the' inflaton potential cannot 
be pursued without making some hypothesis concerning 
the dispersion relation of the fluctuation modes. 

Second, when tensor and scalar fluctuations obey dif- 
ferent dispersion relations, the consistency condition be- 
tween their spectra is modified, as one might have ex- 
pected. What is less trivial is that, in slow roll inflation 
(e2 7^ 0), even if both types of fluctuations are subject 
to the same dispersion relation, the modifications of the 
spectra do not coincide, since the prefactor S in eq. ([35]) 
depends on the value of the parameter fi, see eqs. ([5^ . 
(|40|. and (|46)) . This remark can be important in scenarii 
where H > A. 
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Appendix A: WKB EVALUATION OF THE FIRST 
CORRECTIONS 

We show that WKB waves can be used to reproduce 
the main features of the analytical solution to ex- 
plain their physical origin, and to draw predictions for 
the corrections induced by an arbitrary dispersion rela- 
tion. However the normalisation of the deviations of the 
spectrum cannot be obtained by this method. 

The basis of the argument is that the exact solution 
of eq. ([5]) can always be rewritten as a combination of 
WKB modes: 

u{x) = C{x) v{x) + V{x) iv{x))* , (Al) 



where 



v{x) 



uj{x) 



^2cj(a 



ij{x')dx 



(A2) 



1 2 



F^iH^x- 



Such a decomposition, however, introduces two unknown 
functions and is thus underconstrained. Following section 
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2.4 in [28| we impose: 

- id^u = uj{x)[C{x) v{x) - V{x) {v{x))*] . (A3) 



The unit wronskians of u and v together with 
exactly yields 



and 



\Cix)\'~\V{x)\' = l, 



2uj 



M then 



(A4) 



(A5) 



(C satisfies a similar equation, with C and V interchanged 
and a - sign in the argument of the exponential.) There- 
fore, if in some interval of x, a = dx^juP' — > 0, C{x) 
and 'Dix) are approximately constant in it. We restrict 
our attention to the dispersion relations such that this 
is realized for a; — > oo. In this case the m-mode is sim- 
ply the asymptotic positive frequency WKB solution, the 
asymptotic values of C(a;) and D^x) areC — 1 and T) — 0. 
We also assume that there is a second interval later on 
where this condition is realized. Then V is constant in 
this interval with a non-zero value V f . This value repre- 
sents the probability amplitude of non-adiabatic transi- 
tions and is computed to first order by setting C = 1 in 



V 



f 



2uj 



(A6) 



This integral can be evaluated by a saddle point method 
when the non-adiabaticity is weak. 

We perform this calculation in the particular case of 
the quartic dispersion relation of equation (|10[) . To 
start the analysis we first drop the — ^ term in U!{x), 
which allows us to integrate all the way from x = +oo 
to X = 0. The path of integration can be deformed 
continuously into iM.^. The integral is dominated by 
the contribution from the saddle point x* = 2iX where 



= 0. Using 



I ^1 + ^dx' = Aarcsinh + : 



this yields 



(A7) 



(A8) 



where the constant d can be shown to be equal to 1 in the 
limit A ^ 1 [2^. In this limit, Vf being exponentially 
suppressed, (jA4[) justifies the approximation C ~ 1, and 
the squared amplitude of the mode for 1 ^ \x\ ^ is 
given by: 

|u*"(2;)|2~ |t;(a;)|'(l + 2e-"^cos2(p(x)+ e'^-^) , 

(A9) 

where f{x) is the phase of i;(a;) accumulated from zero to 
X. The important point for us is that 2?/ is exponentialy 
suppressed and not power law suppressed as A ^ 00. 



We can now qualitatively understand the form of the 
rhs in (|17p by the following argument: neglecting the 
term 0{e^'^'^^), before horizon exit, but before the WKB 
approximation completely ceases to be valid, the ratio of 
the instantaneous values of the power spectrum with and 
without dispersion is, 



^0 



\Mx)\' 

|«o(x)P 

UJo{x) 



2 e"''^ cos 2^3(2;)) 



w(x) 



(l-|-2e"''-^cos2^(a;)) , (AlO) 



where wo(x)2 = 1 - ^. To first order in (F'^ - P'^)/P'^ 



we have: 



ujf{x) 



1 - 



2P2 



For the quartic case, one thus has 



uj{x) 



X 

8A2 



(All) 



(A12) 



Thus, eqs. (jA12p and (jA10[) allow us to explain the origin 
of the features observed in (fT7|) : 

i) The exponentially small corrections come from the 
non adiabatic transitions engendered by the dispersion. 
They are thus the result of a global (cumulative) effect. 

ii) The polynomial corrections in p?)) can be viewed as 
an imprint of the different normalisation between the rel- 
ativistic and dispersive modes, around horizon exit (43 |. 
These corrections are thus local in the sense that they in- 
volve only the late time behavior of the modes. The fact 
that these corrections scale as A^^ directly comes from 
the power of A in the non relativistic term in the modified 
dispersion relation. (In addition, since the normalisation 
affects equally the positive and negative frequency WKB 
modes, this explains why the polynomial corrections are 
in factor of both the adiabatic and non-adiabatic terms 
in (Ha).) 

We see in (jA12[) that the normalisation of the super- 
luminous dispersive mode is always smaller than the rel- 
ativistic normalisation, since uj2{x) > luq{x) for all x. 
From this it is tempting to deduce that the deviation 
of the power spectrum is always negative for superlumi- 
nous dispersion relations, and positive for subluminous 
ones. However this conjecture is not correct: the numer- 
ical analysis shows that the sign of the deviations flips at 
a = 3. Therefore the above treatment of the corrections 
to the WKB approximation is only indicative when there 
is mode amplification, that is, strong departure from adi- 
abaticity. Thus it cannot be used to explain why the 
leading signatures of sub- and superluminous dispersion 
found in Figures [7] have exactly the opposite sign. 

However this treatment does contain correct elements. 
To show this, we consider in the next appendix a toy 
model in which the UV evolution is identical to that con- 
sidered here but in which there is no horizon exit, and 
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therefore no mode amplification. From the exact solu- 
tions, we shall see that only the global exponentially sup- 
pressed corrections remain, thereby demonstrating that 
the polynomial corrections arise near horizon exit and 
are linked with the mode amplification (the instability) 
occuring after horizon exit. 



Appendix B: GLOBAL CORRECTIONS AND 
BACK-SCATTERING 



In this appendix, we consider a toy model where the 
evolution is governed by the same mode equation as in 
(fTTjl but with the sign of the non-conformal term ^ re- 



versed: 



di 



1 



X 

4A2 



= 0. 



(Bl) 



we have also slightly generalized the situation by intro- 
ducing the extra parameter fi. 

In the present case there is no mode amplification since 
the square frequency remains positive. We thus compare 
the pair creation probabilities with and without disper- 
sion, and show that only non-adiabatic, exponentially 
suppressed corrections appear. 

It is worth noticing the relationship between equation 
(jBl|) and the wave equation for a massless scalar field in 
a black hole spacetime. When studying the radial wave 
function of s-waves in the momentum representation, af- 
ter having factorized out the dependence on the Killing 
frequency (see for instance the appendix of [2§| for de- 
tails) one obtains in the near horizon region the same 
equation as that of eq. (|Bip . Thus we do not exclude that 
the following exercise be relevant for the trans-Planckian 
signatures in Hawking radiation. 



where Jj-y is the Bessel function of the first kind. 

We want to compute the Bogoljubov coefficients de- 
fined by: 



(B5) 



The squared modulus of /3 gives the particle creation 
probability. Using the identity: 



1 



we get: 



7 sinh^ sinh^ 



sje^^-i - 1 ■ 



(B6) 



(B7) 



2. Particle creation rate in the presence of 
dispersion 



We now consider the wave equation (|B1|1 . The positive 
frequency m-mode is: 



(-4.1 ■ 



where 7 is defined as in the previous subsection and where 
the Whittaker function W is the same as in eq. (IT^ since 
both encode the Bunch-Davies vacuum (the value of the 
second index is different because we have flipped the sign 
of the 1/x^ term in the mode equation). The positive 
frequency owt-mode is: 



^2+ 



A / .x' 



(B9) 



1. Particle creation rate without dispersion 

In the absence of dispersion, the wave equation (|B1 
reads: 



+ i + ^ U = o. 

X 



(B2) 



where we restrain ourselves to /i > 1 / \f% not to have 
mode amplification when a; ^ 0. 

The asymptotic m-mode with positive frequency mode 
with respect to the time 77 oc —x^ for x — > cxd is: 



(B3) 



where is the Hankel function of the first kind, and 



7 = i-\/8/z^ — 1. When a; ^ 0+, the positive frequency 
mode is: 



2 sinh7r7 



\/xJi^(x) , 



(B4) 



where M, a -i is another Whittaker function defined in 

In this case as well, the Bogoljubov coefficients can be 
directly read from the identity (equation 13.1.34 in '2(?|): 



t 2 2 



Hence we get 



r(»7) 

V(i+zl-z|)' 



r(i7) 



(m,.,.)* (BIO) 



The squared modulus simplifies 



1/3: 



1 



-(l + e-'^^e'"'') 



(Bll) 



e27r7 _ 1 

|/3oP X (l + e-"^e"^). (B12) 
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From the second equation, one immediately sees that in 
the hmit A ^ oo the pair creation probabihty in the ab- 
sence of dispersion is recovered. And, as announced, in 
the case when there is no mode amphfication, only ex- 



ponentially suppressed corrections are present, and these 
correspond to non-adiabatic effects. A simplified version 
of this result can be found in a black hole context in 
subsection 5.2 of [13] ■ 
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